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Abstract 

It has been suggested that the modified Steiner action functional has desirable properties for a 
random surface action. In this paper we investigate the scaling of the string tension and massgap in 
a variant of this action on dynamically triangulated random surfaces and compare the results with 
the gaussian plus extrinsic curvature actions that have been used previously. 
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1 Introduction 



The issue of whether a non-trivial continuum Umit exists for gaussian plus extrinsic curvature (GPEC) 
lattice actions of the form 

^= ^(X, -X,)2 + A ^ (l-ri,-n,) (1) 

<ij> Ai,Aj 

on dynamically triangulated random surfaces, is of interest for the construction of well-defined lattice 
versions of string theory ^, ^ as well as for constructing models of membranes in biophysics and 

chemistry. The second term in equ.(l), where the fii are the unit normals on neighbouring triangles, 
is a discretization of the extrinsic curvature and acts as a "stiffness" term. If this term is absent one 
has a gaussian discretization of the basic Polyakov action |^ which gives rise to pathologically crumpled 
surfaces due to the failure of the string tension to scale |^ . The dynamical nature of the triangulation is 
manifested as a sum over triangulations, X^T' canonical (fixed number of points) partition function 

N 

Zn{\) =Y. X{dX,5{Y,X,)eM-S). (2) 

T •' i=l i 

where the delta function is inserted to kill the translational zero-mode, and N is the number of points. 
This means that we have in effect a fluid surface. The GPEC model of equs.(l,2) apparently has a low 
A crumpled phase and a large A smooth phase similar to those displayed by identical models on fixed 
triangulation surfaces ^ where the sum over triangulations in equ.(2) is dropped. The initial work in j^] 
found a pseudo-second order transition on small lattices, but later work |[ || with larger lattices and 
better statistics suggested rather that the transition was higher order, or a crossover phenomenon ^. 

The strongest evidence so far that there is, indeed, a transition comes from the measurements of 
the scaling of the string tension and mass gap carried out in [|[ . An earlier measurement of the string 
tension also found results that were consistent with scaling, but in this the points on the boundary, 
which constituted a large proportion of the total number of points, were physically pinned down j^] 
and the vanishing of the lattice string tension at the critical point was assumed. Although analytical 
calculations suggested that the extrinsic curvature coupling A in equ.(I) was asymptotically free and 
hence that there was no non-trivial theory for finite A in the lattice action, the measurements in were 
strongly indicative of scaling and hence a finite string tension in physical units. This implies a non-trivial 
continuum limit at finite A. 

Pending clarification of the behaviour of the GPEC action on dynamical triangulations, we thought it 
a worthwhile exercise to investigate possible alternative lattice random surface actions in order to see if 
their behaviour was more (or less!) clear-cut. We have already conducted some preliminary simulations 
PH , [l2| of actions containing terms of the form suggested by Savvidy et.al. ||l^, Q that incorporate the 
modified Steiner functional The basic "Steiner" action is just 

Sstemer = - ^ \Xi - Xj\9{aij), (3) 
<i]> 

where 0(aij) = |7r — aij \ and ay is the angle between the embedded neighbouring triangles with common 
link < ij >. This is essentially a coarse discretization of the absolute value of the trace of the second 
fundamental form of the surface, rather than its square which appears in the GPEC action. It was 
observed in that an action containing only this term ran into problems with the entropy of vertices 
in smooth configurations and failed to give a well-defined grand canonical (varying number of vertices) 
partition function. It is a relatively simple matter however to concoct variations on this theme that 
constrain the errant planar vertices somewhat such as 

^1 = ^ E - I + ? E - ^M^^v) (4) 



or even 



2 ^ ■' 2 

<ij> <ij> 



A 



52==El^l + 2 ^ |X,-X,|0(a,,), (5) 

A <ij> 
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where the | A| is just the area of triangle A as seen in the space in which the surface is embedded. In [[l4[ 
another ahernative was suggested in which was modified to 9{aij) = |7r — a.yl^ with ^ < 1 which also 
appeared to improve the convergence of the grand-canonical partition function. 

In [Tl|| we carried out simulations of Si , S2 along with a further permutation combining a gaussian 
term with the Steiner part 

^3 = ^ E - ^^)' + ^ E - ^M^v) (6) 

<ij> <ij> 

and found rather similar behaviour to that seen for the GPEC action on small (72 and 144 nodes) meshes 
- namely peaks in the specific heats for the respective actions. For Si we have 

C = — (< S'ltgjng^ > - < Ssteiner >^) • (7) 

We also see, by visual inspection of the surfaces, a low A crumpled phase and a large A smooth phase. 
Although the gyration radius, a measure of the size of the embedded surfaces, 

^2 = ^^^,^E(^-^.)%^'?. (8) 

where the qi are the number of neighbours of point i, was not monotone increasing with A as for the 
GPEC action this could be explained by noting that the Steiner term, unlike the extrinsic curvature, was 
dimensionful. 



2 Expected Scaling Properties 

Our simulations described above were carried out on boundaryless surfaces with spherical topology. The 
simulations of [Q, which were the most convincing demonstration to date of a non-trivial continuum limit, 
required surfaces with boundaries in order to carry out the scaling measurements of the string tension 
and mass-gap. These are constructed in an ingenious manner using twisted boundary conditions on a 
torus, which we now outline, in order to avoid pinning down a disproportionately large amount of lattice 
at the boundary loops or points. It was observed in |^ that on a torus the sum of vectors Xij along the 
edges of the triangulation on a closed path could take the values 

E{ni,n2) = niEi + n2E2 (9) 

where the vectors i?i,i?2 are constant and the integers ni,n2 denote how many times the path winds 
round the two respective periods of the torus. For non-zero Ei this means that 

X,{ki,k2)^X, + kiEi+k2E2 (10) 

where the ki labelled the particular "copy" of the surface at a given point. The partition function 
in equ.(2) is now dependent on the choice of Ei, Zpf{X) Zn{X, Ei, £2)- Non-zero values of Ei,E2 
correspond to simulating the surface on a frame Ei x i?2. The important point to note is that it is not 
necessary to designate any of the points as boundary points in this procedure. It is thus possible to avoid 
potential poblems with too many points on the boundary for small surfaces. 

A canonical string tension cr(A, N, ) for the system described above is defined by taking Ei — {y, 0, 0) , 
E2 = (0,y,0), Fjv(A,2/2) = -logZjv(A,2/2) and 

-(A,iV,.^) = ^^i^ (11) 

where the translational invariance of Zn means it depends on only the projected area y^. Similarly a 
canonical massgap is defined by choosing Ei = {y, 0, 0), E2 = (0, 0, 0) and 

m{X,N,y)^ . 12) 

oy 
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It is expected that the iV and dependence in the string tension appears as the ratio r — y'^ /N and in 
the massgap as the ratio t = y/N. 

It is actually more natural to define the physical string tension in a grand canonical ensemble (with 
a varying number of points) [p^t , which can be done by taking the Legendre transform of -Fjv(A,y^) 

G{^l,x,y^)^N^, + FN{x,y^) (13) 

where /i is the cosmological constant. For large y^ one expects G(/i, A, y'^) ~ f7(A, with 

a(A,M) = ^-cT(A,r). (14) 

The grand canonical ct(A, /i) is expected to scale as 

a{X,^i)^aa{\)+d{X)^Il' (15) 

where the exponent v governs the scaling of the physical area Aphys — fJ-jiy"^, with /i_R = (/x — ^crit)- It 
is then possible to deduce the expected scaling of the canonical a{X, r): 

a(A,r) ^ao(A) + (Ti(A)r2''/(i-2-). (16) 

The physical string tension (Tphys — (j{X, will be infinite unless ao(A) ~ (A — Acrit)" ^ approach 

a critical point at some Xcrit and this ctq is accessible in a canonical simulation. 
It is also possible to play a similar game with the massgap, defining 

G{^l,X,y) = N^i + FN{X,y) (17) 

which is expected to behave as G'(/i, A, y) ~ to(A, fi)y for large y. In this case we have to(A, /i) ~ /ij^ and 

m(A,i) -£'(A)t"/(i-") (18) 

which is again accessible to a canonical simulation. 

3 Numerical Simulations and Results 

In this paper we apply the methods of |^] to analyse the scaling of the massgap and string tension for one 
of the variant Steiner actions, 5*1. This was chosen because the two terms in the action ^ J2<ij> l^t^-^jl 
and X]<ij> ~ have the same scaling dimensions which simplifies somewhat the choice of 

observables. If we rescale the coordinates Xi — > yX^ in Si we find Si{Xi, y) ySi{X[, 1), which means 
that with the appropriate boundary conditions for the string tension measurements 

dFN{X,r) <Si>-i{N-l) 
= -^y^ = V ^ ^ 

for a surface with N points embedded in 3 dimensions. If we use the boundary conditions that are 
appropriate for the massgap measurements we find 

dFN{X,t) <Si> -3{N - 1) 

to(A, r) ^ = . (20) 

dy y 

We thus simply measure the expectation value of the action with the appropriate choice of frame in order 
to access information about the string tension and massgap scaling. 

In addition, we measure the specific heat as defined in equ.(0) and histogram the output data at the 
various A simulated in order to use the multi- histogram method of Ferrenberg and Swendsen fl^ , which 
allows one to estimate the density of states and hence the specific heat for arbitrary A. We also measure 
the gyration radius, as defined in equ.(^), but a certain amoimt of care is needed with this because of the 
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twisted boundary conditions. We choose to measure the X2 using only the component transverse to the 
frame in the string tension measurements to avoid confusion, and the two components transverse to the 
hne separating the pinned points in the case of the massgap measurements. The autocorrelation times 
for the various observables are calculated in order to ensure that we have reasonable statistics. We also 
measure the various acceptances for the lattice and X moves to check that the Monte-Carlo algorithm, 
which we now describe, is behaving reasonably. 

In order to achieve a reasonable amount of vectorization in the code 64 systems were simulated in 
parallel, with measurements being taken after a sufficient number of sweeps were made to allow them 
to decorrelate. It proved to be convenient to store the link variables Xij rather than the site variables 
Xi which allows the incorporation of the twisted boundary conditions as Xij — Xi ~ Xj + Eij , where 
Eij = n\jEi + nfjE2. The integers riij are non-zero when the link < ij > passes from one of the 
elementary cells in the parameter space (a plane for the torus) to another. Rounding errors during the 
simulation can be kept under control by using the transformations 

Xi Xi + 1]Ei + l^E2 

E,j E^j+llEi+llE2-l]Ei~l]E2 (21) 

where the Vs are arbitrary integers to keep the E'^'s from straying. 

The sum over lattices is effected by carrying out local flip moves on adjacent triangles, forbidding 
flips that lead to degenerate triangulations with less than 3 neighbours per point or with 2-loops. With 
non-trivial boundary conditions the Eij for affected edges must be changed, whereas the Xi are left 
untouched. The coordinates Xi are updated with a simple Metropolis scheme, which does not affect the 
non-trivial boundary conditions. In this paper we report on simulations carried out on relatively small 
surfaces of size 64 and 144 nodes. We have not proceeded to larger surfaces in the current batch of 
simulations because there is a hidden penalty built into the direct transcription of the Steiner action we 
have used in 5i, compared with the GPEC action. Namely, the calculation of 0{aij) requires an inverse 
trigonometric operation, rather than the simple multiplications involved in calculating n-n'm the GPEC 
action. It might be possible to avoid this in further simulations by using some trigonometric function 
with the requisite properties for 9 {9{2tt — a) = 9{a), 9{tt) — 0, 9{a) > ), but this begs the question of 
universality. 

If we now move on to discuss the measurements made for the string tension and massgap with the 
choice of frames Ei — (y,0, 0), E2 — (0, s/,0) and Ei — (y, 0,0), E2 ~ (0,0,0) respectively, comparison 
of Fig.l for the string tension and Fig. 2 for the massgap with Figs. 6, 7 in the first of j|] reveal striking 
similarities. Looking at Fig.l for the string tension first we see that the data points, as expected, fall on 
universal curves as a function of r for a given A until finite size effects set in at small r. Lines are drawn 
to guide the eye through the points coming from the iV = 64 (8 x 8) surface. For large r, just as for 
the GPEC action, we would expect a A independent limit which is what is observed. For small r the A 
dependence becomes more marked and as the A = 3,4 values straddle zero as r — > we can infer there 
is some A value intermediate between these, say Ac where a{X,r) ^ as A — > Ac and r ^ 0. This is 
one of the prerequisites for a finite physical string tension. For A = 4 we see negative <t(A, r) at small r, 
which is due to the repulsion of the vertices, and the value where it is zero corresponds to the equilibrium 
configuration. This is again similar to the behaviour observed in the GPEC action, as are the very long 
autocorrelation times observed in this phase. 

We would expect the results for m(A, t) in Fig. 2 to fall on universal curves for different A with t = y/N 
until finite size effects set in at small t, and this is, indeed, what is seen. We have again drawn lines 
through the = 64 (8 x 8) points to guide the eye. Only the crudest of fits to the scaling law in equ.(18) 
(m ~ D{X)t'^'-^\ where /3(A) = i^(A)/(l - i^{X)) ) is possible with our data, and we find v to be in the 
region of 0.3 for A = 3 and 0.2 for A = 4. This gives an estimate of 1/ at Ac which is lower than that for 
the GPEC action, 0.38 < i' < 0.42. This difference could be due to finite size effects. It has recently been 
pointed out that these may be so strong that they completely mask the presence of the tachyon, and 
hence the real continuum physics, for multiple spin models on dynamical lattices up to very large sizes 
p9[ | , though the effects may well be less severe here as they depend sensitively on the extrinsic Hausdorff 
dimension (of the order of 4 in this case) and probably non-universal constants. We note that our value 
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is consistent with an earlier estimate for the GPEC action in the second of |g| , which was obtained using 
system sizes similar to those studied here. A fit to the scaling law in equ.f? for a{X,r) gives a lower 
estimate for 1/ < 0.2 for both A = 3,4 but the quality of the fits here is even poorer, so the discrepancy 
is not too disconcerting. 

The results for the specific heat are shown in Fig. 3 for some selected r values, where we have not plotted 
the Ferrenberg-Swendsen interpolation between the measured points for simplicity. The Ferrenberg- 
Swendsen interpolation gives a value for Ac that is close to 4. The quality of this interpolation deteriorates 
as the relative size of the frame (ie r or t) is increased, suggesting that finite size effects may be quite 
important in such cases. Finally the time evolution of X2 for r = 1 is shown in Fig. 4 for one of the square 
frames used in the string tension measurements, and is at first sight rather surprising as it decreases with 
time after a disordered start for large A. However, it should be remembered that we measure transverse 
fluctuations to the frame, so if the surface becomes more rigid with increasing A this is what we would 
expect to see as the transverse fluctuations are suppressed. When the frame size is reduced sufficiently 
the customary behaviour of a larger X2 in the large A phase reasserts itself. 

To summarize the numerical results of this paper: for the action 5*1, containing an edge length and a 
Steiner term, the scaling behaviour of both the string tension and the massgap are very similar to those 
seen in the GPEC action. There is some evidence for universal behaviour in the value of the exponent v, 
but our measurements of this are very preliminary. The behaviour of the specific heat peak is consistent 
with that seen in our earlier, smaller scale, simulations of 5*1. Finally X2, both with and without framing 
the surfaces, gives every indication of a smooth or rigid phase at larger A. There is no sign, however, 
of sharper scaling behaviour than is seen with GPEC actions, for Si at any rate. For future work it is 
possible that a subdivision invariant action such as S2 might offer a faster approach to the continuum 
limit [^ . A more judicious choice of 9 from the numerical point of view for any of the Steiner actions 
81,82, 83 might also offer the possibility of more efficient simulations. Nonetheless, the current batch of 
simulations has demonstrated that there is evidence for scaling and hence a non-trivial continuum theory 
for a particular Steiner action, just as with the GPEC action. 
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Figure Captions 

Fig. 1. The canonical string tension a{X,r) is plotted for the various A values to show the scaling with 
r = y'^/N. 

Fig. 2. The canonical massgap m{X,t) is plotted for the various A to show the scaling with t = y/N. 

Fig. 3. The specific heat vs A for some selected r values. Our original small scale simulations of 5*1 [0 
are also included for comparison. 

Fig. 4. The decrease of the perpendicular X2 with time at large A is evident from these measurements 
taken with r = 1.0 
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